Guided plasmons in graphene p-n junctions 
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Spatial separation of electrons and holes in graphene gives rise to existence of plasmon waves 
confined to the boundary region. Theory of such guided plasmon modes within hydrodynamics of 
electron-hole liquid is developed. For plasmon wavelengths smaller than the size of charged domains 
plasmon dispersion is found to be cj oc q 1 ^ 4 . Frequency, velocity and direction of propagation of 
guided plasmon modes can be easily controlled by external electric field. In the presence of magnetic 
field spectrum of additional gapless magnetoplasmon excitations is obtained. Our findings indicate 
that graphene is a promising material for nanoplasmonics. 
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Introduction. Breakthrough progress in synthesis and 
characterization has made graphene [2j a promising ob- 
ject for nanoelectronics. Operation of graphene-based 
transistors [1] and other components would rely on the 
properties of its single-particle excitations - electrons and 
holes. However, one can also envisage a completely dif- 
ferent set of applications which employ collective excita- 
tions, such as plasmons. Currently, plasmon excitations 
in metallic structures are a subject of nanoplasmonics, a 
new field which has emerged at the confluence of optics 
and condensed matter physics with one of the aims be- 
ing the developing of plasmon-enhanced high resolution 
near- field imaging methods [3, HI]. Another objective is 
possible utilization of plasmons in integrated optical cir- 
cuits. However, perspectives of graphene for nanoplas- 
monics are largely unexplored since plasmon modes of 
graphene flakes have not been addressed so far. As our 
results indicate a great amount of control over graphene 
plasmon properties makes it a very promising material 
for applications. 

Fundamentally, the spectrum of collective charge oscil- 
lations reflects the long-range nature of Coulomb interac- 
tion. In conventional two dimensional systems, such as 
those created in semiconducting heterostructures, plas- 
mons are gapless, ui 2 (q) = 2ite 2 nq/m* , with n and m* 
being electron density and effective mass, respectively 
Q. Such oscillations can be treated hydrodynamically. 
In clean graphene at zero temperature the plasmon fre- 
quency, tu 2 cx \Ep\, vanishes with decreasing the doping 
level Ep. It has been argued Q that the interaction be- 
tween electrons and holes in the final state can modify the 
response functions of Dirac fermions and open up a pos- 
sibility for the propagation of charge oscillations at low 
frequencies uj < qv, where v is electron velocity. Still, hy- 
drodynamic (w > qv) analog of conventional plasmons re- 
mains absent unless either temperature is non-zero Q or 
graphene is driven away from the charge neutrality point 
by doping or gating 9]. Expectedly, in both cases plas- 
mon spectrum has the conventional form, w{q) cx q 1 ! 2 . 

In the present paper we investigate spectra of hydro- 
dynamic plasmons in spatially inhomogeneous graphene 



flakes. Realistic graphene samples are typically subject 
to disorder potential and mechanical strain [Toj that lead 
to the formation of charged electron and hole puddles 11 1 
with boundaries between n and p regions being the lines 
of zero chemical potential. Moreover, controlled p-n junc- 
tions can be made with the help of metallic gates 12 1. 
Also p-n junctions can be created by applying electric 
field within the plane of a graphene flake, see Fig. 
The field separates electrons and holes spatially in a way 
that allows control of both the amount of induced charge 
(and thus plasmon frequency) and spatial orientation of 
the junction (the direction of plasmon propagation). 




FIG. 1: Two types of graphene p-n junctions: a) field-induced, 
b) gate-induced. Dot-dashed line indicates boundary between 
electron and hole regions and, correspondingly, the direction 
of plasmon propagation. In case of field-induced junction it 
is controlled by the direction of external electric field Eo- 

Below, we demonstrate that such p-n junctions can 
guide plasmons. We show the existence of charge oscil- 
lations which are localized at the junction and have the 
amplitude decaying with the distance to the junction. 
For wavelengths shorter than the width of the charged 
domains, we find the plasmon spectrum of the form, 



e v Jq\Po\ 



(1) 



where p' is the gradient of equilibrium charge density 
at the junction, v is electron velocity, and n = 0, 1, 2, ... 
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enumerates the solutions. The lowest mode has ao — 
4V27rr(3/4)/r(l/4) w 3.39. 

Below we derive this result and discuss plasmon prop- 
erties for the two types of p-n junctions: electric field 
controlled and gate controlled, as shown in Fig. 1. 

Hydrodynamics of charge density oscillations. We uti- 
lize the hydrodynamic approach to describe the motion 
of charged Dirac fermions. The rate of change of electric 
current density J due to dynamic electric field E follows 
from the usual intra-band Drude conductivity with the 
corresponding density of states [13j, 



j(r,i) 



— a|A*(r)|E(r,t), 



(2) 



determined by the local value of chemical potential /z(r) 
as measured from the Dirac point (positive for electrons 
and negative for holes) . Electric current is related to the 
variation of charge density Sp by means of the continuity 
equation, 



6p{r,t) + V- J(r,i) = 0. 



(3) 



Finally, the variation of charge density produces electric 
field according to the Coulomb law |14| . 



E(r,i) 



r — r 



(4) 



Equations ©-(HI gi ye a closed system for plasmon exci- 
tations in graphene flakes. We apply it to a p-n junction 
created in a strip infinite along the y-axis (direction of 
plasmon propagation) . Using the Fourier representation, 
5p(r,t) = 5p(x)exp(iqy — iuit), and eliminating E and 
J we arrive at the equation for the oscillating part of 
electron density, 



2r , , 2e 2 v d 



po(x)\ d ^ 



dx 



\Po{x)\ 



x / dx'Sp(x')K (\q\\x- x'\) = 0, (5) 

J-d 

Here K is the modified Bessel function and 2c? is 
the width of graphene flake. Within the Thomas- 
Fermi approximation equilibrium charge density po(x) 
is related to the chemical potential via po{x) = 
—sgn(p)ep 2 (x)/Trh 2 v 2 (electron charge is taken to be 
— e). This follows from the condition that the electro- 
chemical potential p(x) — e<j){x) is constant throughout 
the system. The solutions of Eq. ([5]) will now be consid- 
ered for large and small plasmon momenta separately. 

Short wavelength, q 1/d. In this case the decay 
of plasmon density Sp(x) occurs over a distance much 
smaller than the width of the system and the limits 
of integration in Eq. (0 can be extended to infinity. 
Assuming (cf. Eq. (fTT]) below) the linear dependence, 
Po(x) — p' x, we observe that the integro-differential 



equation §5§ acquires obvious scaling property. Intro- 
ducing the variable £ = qx we arrive at the plasmon 
spectrum in the form (p}, with dimensionless constants 
a n determined from the eigenvalue problem: 

x deV n) (0«o(ie-c'i) = on<y/» (n) (0- (6) 



Interestingly, this integro-differential equation allows a 
complete analytic solution, though the detailed analysis 
is beyond the scope of this paper. Our main findings 
are as follows. Solutions are enumerated by n = 0,1,2... 
with even/odd numbers corresponding to even/odd den- 
sity profile, M n) (-0 = (-l) n fy (T,) (£)- Surprisingly, 
eigenvalues are doubly-degenerate and given by 



2n + l 3-7--(4n-l) 
a2 "- ao 4^TT'l.5-.(4n~3)' a2n+1 



«2r, 



(7) 



At large distances all modes have exponential depen- 
dence, <5p (n) (0 - e-l ? l, while at |£| < 1 even and 
odd solutions exhibit different behavior, 5p( evcn ^ ~ 1 — 
const and (5p( odd ) - sign($)/ VTCf- The fi rst pair 
of solutions (belonging to the lowest eigenvalue ao) in 
the Fourier representation Sp^(k) = J d^Sp^ (£)e lfe ^ 
acquires a simple form: 



Sp {0) (k) oc 



(1 + fc 2 ) 3 / 4 ' 
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(1 + fc 2 ) 3 / 4 ' 
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Long wavelength, q <C 1/d. In contrast to the above 
result ([lj plasmon spectrum at small q is sensitive to a 
specific realization of the p-n junction. We address the 
long- wavelength behavior of plasmons in field controlled 
junctions. We expect this case to be of more interest, 
in addition it allows a more complete description. Be- 
fore analyzing plasmons in this structure, we discuss the 
equilibrium density profile. As shown in Fig. la the flake 
of width 2d is placed in external electric field Eq applied 
along the x-direction. The equilibrium density distribu- 
tion p(x) is found from, 



x ~\~ x 1 

dx'p (x') In = 0, 

\x — x'\ 



, flV IT . ... 

E x+sgn{x) — x -\p (x)\+2 

.In !■' ■' 

(9) 

where it is used that po(x) = —po(—x). Prior to solv- 
ing Eq. © it is instructive to analyze validity of the 
semiclassical approach. The first condition implies that 
the change of the electron wavelength is smooth on the 
scale of itself, d/dx(hv/p) <C 1. Estimating fi(x) ~ eE x 
we obtain that the distance to the p-n junction line 
(x = 0) should exceed the characteristic electric field 
length Ie — y s/Eq <c x. The second condition requires 
that the electron wavelength is small compared with the 
width of the system, d 3> Tiv / p. Noting that in graphene 
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hv ~ e 2 we can rewrite this second condition simply as 
Ie *C d. Thus, the Thomas- Fermi equation ([9]) for the 
equilibrium charge density and the hydrodynamic equa- 
tion ([5]) for its variation are applicable as long as 



l E <^d, q<l/h 



(10) 



However, the ratio of q and 1/d can be arbitrary. For a 
moderate external electric field ~ 10 4 V/m the value of 
electric length Ie ~ OApm and the first of the conditions 
(fl"0| is satisfied easily for micron-sized samples. 

Analytic solution of Eq. (j9|) is possible when the second 
term is small, in which case the charge density is 



EnX 



Vd 2 



(11) 



Substituting this expression back into Eq. ^ we ob- 
serve that the second term is indeed negligible as long 
as x 3> 1%/d. This is assured whenever the condi- 
tions (fTU|) are satisfied. It is also worth pointing out 
that Eq. (|11[) justifies the linear approximation for the 
charge density used in deriving Eq. ([1]) for <; > 1/d, with 
p' /e = l/(l%d). 

We now turn to the analysis of plasma oscillations 
propagating on top of the density distribution, Eq. (fTTJ) . 
For small plasmon momenta, q -C 1/d, electric field ex- 
tends beyond the width of the flake and the equation ([5]) 
needs to be supplemented with the boundary condition, 
which ensures that electric field (and thus the current) 
vanishes at the edges, x — ±d: 



dx 



x ± d 



0. 



(12) 



The spectrum of the lowest symmetric mode can be most 
easily found by integrating Eq. ([5]) across the width of 
the flake. The first term in the brackets will then van- 
ish exactly due to the boundary condition ([12]) , The 
remaining integral can now be calculated to the log- 
arithmic accuracy with the help of the approximation 
Kf)(q\x — x'\) — — \aq\x — x'\: 



dx 



Mf, , , ,.v 2dr 2 (3/4) 

In (q\x - x ) « — =— In (qd). 

e IevK 



(13) 

Eqs. ([5]) and (ITBl combine to give the equation, [uj 2 — 

^oil)] f-dd x 3p( x ) = 0> that yields the dispersion of the 
gapless symmetric plasmon, 



W g( ? )=I*(3/4)^ g a ln(W), 
irnlE 



(14) 



reminiscent of the plasmon spectrum in quasi-one- 
dimensional wires, The remaining modes, n > 1, are 
gapped. For these modes J^ d dx5p(x) = and simple 
procedure of integrating Eq. (JSJ over the width of the 



flake is not useful. Instead, the equation for the n-th fre- 
quency gap can be obtained by setting q = in Eq. ([5]). 
We observe that 



Wn(0)=A 



e v 



(15) 



' hlEd' 

where f3 n are the eigenvalues of the equation, 

2 * vw\_f^_^m = ^ M{() iw) 



^dt(i-e) i/4 



The zeroth mode /3o = 0, see Eq. (Tl4|) . is found ana- 
lytically: Sp^ oc 1/y/l — C 2 - It describes charge dis- 
tribution in the strip in response to a (uniform along 
x direction and smooth along y-direction) change of its 
chemical potential [16j . Other solutions of Eq. ()16|) are 
found numerically: 



181 = 1.41, /3 2 =6.49, /3 3 =6.75,.. 



(17) 



With increasing n the eigenmodes of integro-differential 
equation (fT6|) oscillate faster, but in general do not follow 
the oscillation theorem familiar from quantum mechan- 
ics. In particular, the solutions with n — and n — 3 are 
even while n = 1, n = 2 are odd [17| . 

Finally, we mention the case of a gate-controlled p-n 
junction, Fig. lb. The equilibrium density profile is linear 
near x = and saturates for large \x\ [18j . Eq. ([1]) is still 
applicable for q > 1/d. In the limit q < 1/d one should 
take into account the screening of long-range Coulomb 
interaction by metallic gates. In this case the logarithm 
in the spectrum of the gapless plasmon disappears, and 
the lowest mode Eq. (fl"4)) becomes sound-like. 

Magnetoplasmons. If external magnetic field B is ap- 
plied perpendicularly to the plane of graphene the plas- 
mon spectra acquire new modes. The equation of motion 
@ should now be modified to include the Lorentz force, 



9 2 
• e ev 

J(r,t) = — aMa:)|E(r,t) 

nh cp(x) 



J x B. 



(18) 



The relative coefficient between electric and magnetic 
terms in this equation follows from the expression for 
the Lorentz force acting on a single particle. The last 
term has opposite sign for electrons and holes. Note that 
the frequency of cyclotron motion lob(x) — ev 2 B/cp{x) 
in graphene p-n junctions is position-dependent. The 
remaining equations ((3])-(j4|) are intact in the presence of 
magnetic field. The boundary condition requires now the 
vanishing of the normal component of electric current at 
the boundary, rather than simply vanishing of the elec- 
tric field, as in Eq. (fl"2"f . Eliminating J and E we arrive 
at the generalization of equation ([5]), 

«*> + v{' 22 -> 2 >'-^} 

x/ dx'6p(x')K (\q\\x-x'\)=0, (19) 

J-d 
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where Z(x) — \/i(x)\/(oj%(x) — w 2 ). 

The most interesting effect described by Eq. (Tl9t is 
the appearance of a set of new modes, chiral magne- 
toplasmons, similar to those considered in Ref. [l9[ for 
conventional 2D electron systems with smooth bound- 
aries. To find their dispersion in strong magnetic fields, 
when u) <C lob(x) (the exact condition is given below), 
one should retain only the second term in Eq. (fT9|) . 
Noticing that (ujbZ)' = jrl B p' (x)/e = nl B /(l E d), where 
Ib = \/Tic/eB is the magnetic length, we arrive at the 
integral equation 



2c q dpp(x) 
B oj dx 



dx'Sp(x')K (\q\\x - x'\) = Sp{x) (20) 



Since Kq is positive, propagation of magnetoplasmons 
with q > is quenched, indicative of their chiral property 
[13]. As seen from Eq. (|20|) . the plasmon density Sp(x) is 
concentrated where p' Q (x) is the strongest. The derivative 
of the charge density in field-induced junctions (fTT|) fea- 
tures strong singularity near the edges of the flake. Thus, 
low-frequency magnetoplasmon spectrum is strongly de- 
pendent on the microscopic regularization of this singu- 
lar behavior and is, therefore, beyond the scope of the 
Thomas-Fermi approximation used throughout this pa- 
per. 

The gate-induced junctions, however, allow a rather 
simple analytical description of these modes if we ap- 
proximate that p' Q (x) — e/l 2 E d for |a;| < d and p' Q (x) = 
for > d. The oscillating density 8p(x) then vanishes 
for |x| > d. The solution inside the strip, \x\ < d, can 
be easily found for q ^> l/d, where one can assume the 
range of integration in Eq. (|20|) to be infinite. The eigen- 
functions of Eq. (|20|) are simply given by sin [q± (x + d)} , 
with the values of q± — ixnj2d determined from the con- 
dition, Sp q (±d) = 0. The spectrum of magnetoplasmons 
is then found to be, 



2 7 re 2 Z 2 , 

w„(«) = - 



M%d 



yjq 2 + 7r 2 n 2 /4d 2 ' 



re=l,2... (21) 



The magnetoplasmon spectrum (|21l) is derived under 
the assumption that magnetic field is strong, uj E (d) 3> w, 
which implies that Ib <SC Ie- In order to neglect the first 
and third terms in the brackets in Eq. (TH)1) one has to 
ensure that q <C (Ie/Ib) 4 /d. This condition might turn 
out to be more or less restrictive than the hydrodynamics 
condition q <C 1/Ie, depending on the particular value of 
the ratio Ib/Ie- Note that the smallness of this ratio is 
not in contradiction to the non-quantized description of 
electron motion in magnetic filed. The latter is valid as 
long as the filling factor is large, eEd ^> uJs{d), which 
means that Ib l E /d. For magnetic field ~ IT, and 
Ib ~ 25nm, using the estimate below Eq. (|10[) that Ie ~ 
400nm we conclude that the width of the flake should 
exceed d > 10pm. The magnetoplasmon modes (j2"T1) are 
~ (Ib/Ie) 2 slower than electrons. Note that these modes 



are undamped since single-particle excitations cannot be 
induced at frequencies below cyclotron frequency u>b- 

Conclusions. Graphene p-n junctions are among the 
most simple and promising applications of this material. 
Single-electron properties of p-n junctions have been ex- 
tensively studied. In the present paper we investigated 
their collective excitations both with and without mag- 
netic field. We anticipate that plasmon modes will be 
crucial for the optical response of graphene nanostruc- 
tures and realistic samples containing electron-hole pud- 
dles. High degree of experimental control should make 
them of special interest to nanoplasmonics and electron- 
ics. Among the most promising applications of plasmons 
in p-n junctions we envisage a possibility of a "plasmon 
transistor" In particular, by simply switching the 
direction of electric field from across the flake to along 
it (and back) the propagation of plasmons can be facil- 
itated (or prevented). In addition, as follows from the 
above Eqs. ([1]), (fTT|) . the plasmon velocity can be con- 
trolled with simple change in the magnitude of electric 
field. This is in a sharp contrast to plasmons in metal- 
lic nanostructures, whose spectra are typically fixed once 
devices are fabricated. 
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